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CO ■ Abstract. We study effects of delay in diffusively coupled logistic maps on the Cayley tree networks. We 

find that smaller coupling values exhibit sensitiveness to value of delay, and lead to different cluster patterns 
of self-organized and driven types. Whereas larger coupling strengths exhibit robustness against change 
in delay values, and lead to stable driven clusters comprising nodes from last generation of the Calaye 
$H ' tree. Furthermore, introduction of delay exhibits suppression as well as enhancement of synchronization 

Oh, depending upon coupling strength values. To the end we discuss the importance of results to understand 

■ conflicts and cooperations observed in family business. 
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1 Introduction larity of Cayley trees makes analytical understanding or 

origin of driven patterns easier to understand using Lya- 
£H Many real- world networks display local co-ordination among p Unov function analysis 

^ , nodes leading to cluster synchronization dJUEE] . Forma- r av i PV Tr( ™ nroviHe a sirrmle model to Ho exact anal 
„ . i • i i i ii i • % ^ayiey Liees pioviue a simple iiiuuci 10 uo exad anai- 

Cj ' tion ol clusters, which are based on the dynamical proper- . c , , c , • j , , rofii <. t j i 

i " i . ,. „ , ' . „ , i ,, i ysis for stability of synchronized states 1201, to study lo- 

" tics of the coupled system, typically depends on the under- ,. ,. .... . p— *~ t A . 

. . , , f r^i • , i i , ,i , cahzation criteria m impurity atom 1211 , to derive expres- 

lying network structure, the interplay between the struc- c ,• j c rp i-u-i-j. nr^n 

CN • r i , ii i S10n i° r magnetization and zero held susceptibility 1221 . 

>ture and the dynamics of complex networks has been a , r>- i • n ■ j. j i <-i i ± a i^ 

. . , J ,. -,,i , f, „ etc.. Biologically oriented work on Cayley tree networks 

locus of intense research interest m last decades 5 . .bur- . , , , ,. e ■ , , , , 

\T .... include modeling of immune networks with antibody dy- 

\r~\ • thermore, delay naturally arises m extended systems due . r?™ T , ^ , , , , , 

*■ ' , , „ !, ' , r . r J . . . im t-, namics 23 . In a recent paper, Cayley trees have been used 

CO to the hmte speed of information transmission 6 . tor ex- , . , ^. , . , , . r?r-n 

1_ ■ ^ , ™ ... to investigate hsosc-Fmstem condensation 24 . 

~ ) ■ ample, m neural networks, propagation delays ol electrical 

£f*l ' signals connecting different neurons and local neurovascu- 
, lar couplings lead to time delays • A delay may give 
CO ' rise to many new phenomena in dynamical systems such as 
| oscillation death, enhancement or suppression of synchro- 
^ ■ nization, chimera state, etc [Ul fTUlfTTlfT^lfTB] . The existence 

of delay can completely change the behavior of a system as We use well known de i ayed C0U pled maps model 
observed for undelayed case [5] ■ What follows that time de- 
lay might be deliberately implemented in order to achieve 
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desired functions such as secure communication Q3I and /. . -,\ n \r f /,u £ \ • , , ,, , -, . 

i ii-i ' r ' i x i(t + 1) = (1 ~ £ ) J \ x i\t)) H / AijQ\Xj(t — r)J (1) 

to control neural disturbances, e.g., suppression of unde- fcj f-^ 

sired synchrony of firing neurons in Parkinson's disease or 

epilepsy nnmriT]. 

Our recent work demonstrated that delay plays a cru- Here h = J2j=i A ij 1S degree, and Xi(t) is the dynamical 

cial role in formation of synchronized clusters and mech- variable of the i — th node (1 < i < N) at time t, A 

anism behind the synchronization. We presented results is the adjacency matrix with elements taking values 

for cluster formation in delayed coupled maps on 1-d lat- 1 and depending upon whether there is a connection 

tice, small-world, scale-free, random and complete bipar- between i and j or not. The delay r is the time it takes 

tite networks [TH] ■ In this paper we investigate delay- induced for the information from a unit to reach its neighbors and 

cluster patterns in diffusively coupled logistic maps on be processed. The function f(x) defines the local nonlinear 

Cayley tree networks. map and the function g(x) defines the nature of coupling 

The Cayley tree is an infinite dimensional regular graph between the nodes. We present the results for the local 

with an idealized hierarchical structure |19j . Its rich hier- dynamics given by the logistic map f(x) — fix(l — x) and 

archal structure turns out to be an ideal model network to for diffusive coupling g(x) — f{x). We take fi = 4, for 

investigate driven patterns in details. Furthermore, regu- which logistic map exhibits chaotic behavior. 
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Fig. 1. (Color online) Phase diagram demonstrating different 
values of (a) finter and (b) f intra in two parameter space of e 
and r ((a) and (b) ) with N = 127, < k >= 2. Local dynam- 
ics is governed by logistic map f(x) = 4a;(l — x) and coupling 
function g(x) — f(x). The figure is obtained by averaging over 
20 random initial conditions. The color-scale encoding repre- 
sents values of /inter and f intra- The regions, which are black 
in both graphs (a) and (b), correspond to states of no clus- 
ter formation. The regions, where both subfigures have gray 
shades (yellow) , correspond to states where clusters with both 
inter- and intra-couplings are formed. The regions in (a), which 
are lighter as compared to the corresponding e and r values 
in (b), refer to dominant D phase synchronized clusters states, 
and the reverse refer to dominant SO phase synchronized clus- 
ters state. White (light yellow) regions in (a) and (b) refer to 
ideal D and ideal SO clusters state respectively. The regions, 
which are dark gray in (a) and black in (b) or vice-versa, corre- 
spond to states where very few nodes are forming the cluster. 

3 Phase synchronization and synchronized 
clusters 

Synchronization of coupled dynamical systems is denned 
by the appearance of some relation between the functional 
of different dynamical variables. The exact synchroniza- 
tion corresponds to the situation where the dynamical 
variables for different nodes have identical values. The 
phase synchronization corresponds to the situation where 
the dynamical variables for different nodes have some def- 
inite relation between their phases [25] . We consider phase 
synchronization as defined in [51126). 

Phase synchronized clusters: Let ni and rij denote the 
number of times the variables Xi(t) and Xj(t),t = 1, 2, . . . T 
for the nodes i and j, show local minima during the time 
interval T. Let denote the number of times these local 
minima match with each other. Phase distance between 
two nodes i and j is given as dij = 1 — 2riij/(ni + rij). 
Clearly, dij — when all minima of variables Xi and Xj 
match with each other dij — 1 when none of the minima 
match. Nodes i and j are phase synchronized if dij = 0. A 
cluster of nodes is phase synchronized if all pairs of nodes 
of the cluster are phase synchronized. 



(a) Cluster node synchronizes with other nodes and forms 
a synchronized cluster. Once this node enters a synchro- 
nized cluster it remains in that cluster afterwards. 

(b) Isolated node does not synchronized with any other 
node and remains isolated all the time. 

(c) Floating node keeps on switching intermittently be- 
tween an independent evolution and a synchronized evo- 
lution attached to a cluster. 

The study of relation between the synchronized clus- 
ters and the coupling between the nodes represented by 
the adjacency matrix exhibits following two different phe- 
nomena of cluster formation. 

(f ) Self-organized clusters: The nodes of a cluster can be 
synchronized because of intra-cluster couplings. We refer 
to this as the self-organized (SO) synchronization and the 
corresponding synchronized clusters as SO clusters. Ideal 
SO synchronization refers to a state when clusters do not 
have any connection outside the cluster, except one. Dom- 
inant SO synchronization corresponds to the state when 
most of the connections lie inside the cluster. 
(2) Driven clusters: The nodes of a cluster can be synchro- 
nized because of inter-cluster couplings. We refer this as 
driven (D) synchronization and the corresponding cluster 
as D cluster. The ideal D synchronization refers to the 
state when clusters do not have any connections within 
them, and all connections are outside. Dominant D syn- 
chronization corresponds to the state when most of the 
connections lie outside the cluster and very few inside. 

To get a clear picture of self-organized and driven be- 
havior we consider two quantities f intra and finter as mea- 
sures for intra-cluster and inter-cluster couplings as fol- 
lows: 

r Nintra ? A ^ inter /~\ 

J intra — ) J inter (2j 

where Ni ntr a and A^„ ter are the numbers of intra- and 
inter- cluster couplings, respectively. In Ni nter , coupling 
between two isolated nodes are not included. 

We define one more important state which forms basic 
backbone of the present investigation. 
Cluster patterns: A cluster pattern refers to a particular 
phase synchronized state, which contains information of 
all the pairs of phase synchronized nodes distributed in 
various clusters. A cluster pattern can be static or dynam- 
ical. Static pattern has all nodes fixed, except few floating 
nodes, in a cluster with respect to change in time, de- 
lay value or initial condition. Dynamical pattern exhibits 
changes with time evolution, or with initial condition or 
with change in delay value. A change in the pattern refers 
to the state when members of a cluster get changed. Fur- 
thermore, patterns can be of D or SO type, which respec- 
tively refers to a particular D or SO phase synchronized 
state. 



4 Mechanisms of cluster formation 

Depending upon the asymptotic dynamical behavior the 
nodes of the network can be divided into the following 
three types [2"T] . 



5 Numerical Results 

We evolve Equation ([I} starting from random initial con- 
ditions, and study the dynamical behavior of nodes after 
an initial transient. We study phase synchronized clusters 
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Fig. 2. (Color online) Schematic network diagrams illustrat- 
ing different cluster patterns observed for different delay val- 
ues at lower coupling strength region. The examples are for 
N = 31, < k >— 2 and e = 0.16. The closed circles of same 
shade (color) imply that the corresponding nodes are phase 
synchronized (i.e. C-y = 1), and the open circles imply that 
the corresponding nodes are not phase synchronized. The D 
chaotic clusters for r = 0, 2 and 4. The SO periodic clusters 
for r = 1, 3 and 5. 



for 100 time steps after the initial transient, and calculate 
values of fmter and fintra as described earlier. 

Phase diagrams Figs.(QJi) and (JTJj) are plotted for net- 
work size N = 127 and average degree < k >— 2. Un- 
delayed Cayley tree exhibits dominant D clusters in the 
range 0.12 < e < 0.19 with a periodic dynamical evolu- 
tion. With a further increase in coupling strength, there 
is no phase synchronization till e = 0.4, after which dom- 
inant D clusters are obtained as elucidated by light gray 
(yellow) regions in Fig. ((1^,) and dark gray (red) regions 
in Fig. ((TJd). At very high coupling strengths, persistence 
of light gray (yellow) regions in (JTJi) and appearance of 
black regions in (QJi) indicate ideal D clusters. 

On introduction of a delay r = 1 in the evolution 
equation Eq. [TJ after very small coupling values for which 
there is no phase synchronization (black color for the sub- 
figures ([Th.) and ([lb)), SO phase synchronized clusters 
are formed for 0.12<e < 0.19 as elucidated by white re- 
gions in ((TJ)). As coupling strength increases, in the range 
0.36 < e < 0.42, where undelayed system exhibits no or 
very less cluster formation, delayed evolution manifests 
dominant D clusters as depicted by gray (yellow) regions 
in Fig. (QJi). With a further increase in coupling strength, 
appearance of gray (yellow) regions in Fig. Pi) and black 
regions in Fig. (JTJd) indicate ideal D clusters. 

For r = 2, appearance of white (light yellow) window 
in Fig. (JTJl) and corresponding window in (QJ>) with dark 
gray (red) to black shades indicate formation of dominant 
and ideal D clusters respectively in lower coupling values. 
This description is same as observed for r = 0. Larger 
coupling strengths lead to ideal D clusters as elucidated 
by black and light gray (yellow) regions in Figs. (HJd) and 
([Tk.) respectively. 

For a further increase in delay, at lower coupling strength 
odd delay values exhibit similar behavior as observed for 
r = 1, while even delay values manifest similar behavior 
as observed for r = and t = 2. At higher coupling val- 
ues, coupled dynamics for all delays demonstrate either 
ideal or dominant D clusters. 



5.1 Delay-induced change in mechanism of cluster 
formation: 

Above discussions indicate that at the lower coupling val- 
ues, change in delay values are related with the change in 
mechanism behind the cluster formation. Odd delay val- 
ues lead to ideal or dominant SO clusters, whereas even 
delay values are associated with ideal or dominant D clus- 
ters. Figs. ([2ji), fl2j:) and (03) illustrate that for r = 0, 
r = 2 and r = 4, nodes in alternate generations are syn- 
chronized with each other, except few cases where nodes 
in two consecutive generations (parents and children) too 
exhibit synchronization. Figs. $jpi), ([Hi) and fl2J) demon- 
strate that for t = 1, t = 3 and r — 5, either one single 
cluster is formed spanning all nodes, or several clusters 
are formed with clusters consisting of nodes in consecu- 
tive generations. 

The cluster patterns observed here are dynamical with 
respect to intimal condition as well as with delay value, 
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Fig. 3. (Color online) Schematic diagrams illustrating delay- 
induced driven patterns for larger coupling values. The exam- 
ples are for N = 40, < k >= 3 and e = 0.37. The closed circles 
of same number (same color) imply that the corresponding 
nodes are phase synchronized (i.e. Aij = 1), and the open 
circles imply that the corresponding nodes are not phase syn- 
chronized. The D chaotic clusters for r = 1, 2 




Fig. 4. (Color online) Schematic diagrams illustrating effects 
on delay on phase synchronized patterns. The examples are 
for N = 30, < k >= 3 and e = 0.7. Closed circles having 
same number (same color) imply that corresponding nodes are 
phase synchronized in same cluster, and open circles imply that 
corresponding nodes are not phase synchronized. 

but for a particular value of delay the phenomenon be- 
hind the synchronization in cluster-pattern is static and 
same parity of delay leads to same phenomenon of cluster 
synchronization. The dynamical evolution in this range of 
coupling strength is periodic for all delay values. 



5.2 Delay-induced driven patterns 

As described in earlier sections, for coupling range 0.35 < 
e<0.42, undelayed coupled maps do not exhibit cluster 



formation, whereas delayed evolution leads to dominant 
D clusters. In order to explain different these dynamical 
cluster patterns clearly we make schematic diagram of dy- 
namical clusters in Fig. ([3]). For undelayed evolution there 
is no co-ordination between any pair of nodes, and hence 
there is no cluster pattern as depicted by all empty circles 
in Fig. <[3^l) . Introduction of delay induces co-ordination 
between nodes in the same sub-family of the last genera- 
tion, as depicted by different clusters in Fig.([3}D). Further 
change in delay value does not have any measure impact 
on synchronized clusters state. These cluster patterns are 
stable with respect to time evolution, initial condition as 
well as change in delay value. 

For a further increase in e, delay destroys the co-ordination 
between the nodes which are connected, giving rise to 
ideal D patterns with again only last generation nodes 
being synchronized in several clusters (Fig. (Hp)). These 
patterns are too stable with respect to change in initial 
condition or change in delay value. Furthermore, delayed 
as well undelayed dynamical evolution in this range are 
associated with the chaotic evolution. 



6 Lyapunov function analysis 

As demonstrated above, while undelayed evolution in mid- 
dle coupling strength exhibits some synchrony between 
children in last generation and their parents, and hence 
giving rise to SO clusters, delayed evolution yields only 
D cluster indicating loss of synchrony between parent and 
children. Introduction of delay destroys the synchroniza- 
tion between parents and children while keeping the co- 
ordination between children unaffected. In order to under- 
stand this behavior let us write down Lyapunov function 
for a pair of synchronized nodes as, 

V ij (t + l) = [(l-e)(f(x i (t))-f(x j (t))) + 

N 

=]v — — A ikg(x k (t - t)) - 

l^k=l Ai i k=l 

N 

^A tfcg (x fc (t-r))] 2 

l^k=l A ij k=l 

Lyapunov function takes following simple form for two 
nodes originated from the same parent node, 

V ij (t + l) = [(l-e)(f(x i (t))-f(x j (t)))} 2 

Above equation does not have any delay term. All this 
leads to the conclusion that if a pair of nodes originated 
from same parent are synchronized, the introduction of 
delay would not affect the synchrony between them. In- 
troduction of delay only affects the connected nodes such 
as parent and children, as it removes the common term in 
their evolution equations, and hence may be a reason be- 
hind the destruction of co-ordination between them |18j . 
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7 Conclusion 

We have studied delay-induced patterns in coupled maps 
on Cayley tree networks. We demonstrate that different 
delay values manifest different cluster patterns at lower 
coupling values, where change in delay not only leads to 
a completely new pattern but also associated with differ- 
ent mechanisms behind the cluster formation. In middle 
coupling range, delayed evolution always exhibits D clus- 
ters. Though we always find either ideal or dominant D 
clusters in this range, the role of delay on evolution and 
co-ordinations among nodes are very different for differ- 
ent coupling values. For some coupling values where un- 
delayed evolution does not exhibit any synchronization, 
a introduction of delay enhances synchronization between 
children in last generation yielding D clusters, whereas for 
some coupling values delay destroys existing synchrony be- 
tween parents and children while keeping children coordi- 
nated, again giving rise to D clusters. These delay induced 
D clusters are stable with respect to the change in delay 
value and consists only last generation nodes. Lyapunov 
function analysis provides some hints about formation of 
stable D clusters in this region. 

The model which we have considered here demonstrates 
that lower coupling strength in general favors synchro- 
nization in various generations in the family, as indicated 
by larger cluster size and by large number of nodes (al- 
most all) participating in clusters. These clusters are sensi- 
tive with respect to external conditions such a time delay. 
Where as higher coupling strength leads to very drastic be- 
havior such as formation of stable D clusters comprising 
only last generations. The origin of these stable D clus- 
ters for last generation nodes can be very well understood 
for Calaye three, where coupling environment of the last 
generation nodes belonging to same sub-family remains 
same, and hence gives rise to a stable driven cluster [18], 
whereas nodes originated from the same parent in any pre- 
vious generation can not have same coupling environment 
unless their all children are synchronized with each other. 

Coupling strengths can be interpreted as closeness or 
bonding among family members, for example lower cou- 
pling strengths can correspond to a situation where mem- 
bers live in nuclear families and do not share much de- 
tails apart that they belong to a same big family. Where 
as larger coupling strengths can be treated as a situation 
where all members of a family live together as a joint fam- 
ily [2 8) . Our results may be used to understand conflicts in 
brothers running a successful family business [29] . which 
on very simple terms can be attributed to the conflicts be- 
tween their children (as shown in Figs. (3) and (4)), where 
as lower coupling strength keeps a warmer relation lead- 
ing to cooperation in family (as seen for Fig. (2)). Lower 
coupling strength here can be considered as separating the 
business of siblings and cousins, which have been proven 
to increase cooperation between them [30] . 

To conclude, we demonstrate that delay in spatially 
extended systems may lead to a completely different re- 
lation between the functional clusters and topology than 
exhibited by undelayed evolution, and hence provide an 
additional step towards ongoing research attributed to 



understand relation between these two. Since delay has 
already been emphasized to be important for many real 
world networks [5J, the results presented in the paper is 
important to understand various different behaviors ex- 
hibited by these systems. Furthermore, observation of dif- 
ferent cluster patterns as a function of delay may shed 
some light in understanding conflicts or cooperation in 
family business [3T] , 
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